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INTRODUCTION 


Tf  target  acceleration  could  be  measured,  it  is  often  thought  that 
it  could  be  advantageously  used  to  guide  a missile.  However,  it  is  not 
clear  as  to  how  to  determine  target  acceleration,  exactly  what  to  do 
with  it  once  one  has  it,  and  what  can  be  gained  by  having  it.  In  a 
limited  sense,  this  report  tries  to  answer  these  questions. 


In  the  recent  past  significant  research-"  -w  has  appeared  in  open 
literature  describing  terminal  guidance  laws  purporting  to  be  improve- 
ments (in  the  sense  of  reducing  terminal  miss  distance)  over  proportional 
navigation.  The  derivation  of  such  guidance  laws  is  often  based  on 
optimal  control  theory  (see  Footnote  4)  and  makes  certain  assumptions 
about  missile/target  geometry  and  velocities  in  order  to  arrive  at  a 
linear  model  more  amenable  to  analysis.  In  many  cases  the  "optimal" 
guidance  law  consists  of  proportional  navigation  plus  a term  proportional 
to  the  target  acceleration,  for  example 


»(“  - k) 


1U-S.  Army  Missile  Command,  Redstone  Arsenal.  Optimal  Controllers  for  Hominy  Missiles,  by  G.  Willems.  Huntsville, 
Alabama,  11  September  1968.  (Redstone  Arsenal  Report  RE-TR-68-15. 

2U.S.  Army  Missile  Command,  Redstone  Arsenal.  Optimal  Controllers  for  Homing  Missiles  with  Two  Time  Constants,  by 
G.  Willems.  Huntsville,  Alabama,  October  1969.  (Redstone  Arsenal  Report  RE-TR-69-20. 

3E.  I.  Axelband  and  F.  W.  Hardy.  "Quasi-Optimum  Proportional  navigation,."  in  Insu.  Klee,  electron.  Eng. . Trane., 
Automat.  Contr.,  Vol.  AC-15,  No.  6 (December  1970). 

4A.  e.  Bryson  and  Y.  Ho.  Applied  Optimal  Control.  Blclsdcli  Publishing  Company,  Waith.iu,  Massachusetts,  1969. 

5C.  J.  Nazarof f . "An  Optimal  Terminal  Guidance  Lav,"  in  Inst.  Klee.  Electron.  Eng.,  Tr\ms.,  Automat.  Contr., 

Vol.  AC-21,  No.  3 (June  1976). 

^J.  L.  Speyer.  "An  Adaptive  Terminal  Guidance  Scheme  Based  on  an  Exponential  Cost  Criterion  with  Application  to 
Homing  Missile  Guidance,"  In  Inst.  Elec.  Electron.  Eng.,  Trans.,  Automat.  Contr.,  Vol.  AC-21,  No.  3 (June  1976). 

7The  Analytic  Sciences  Corporation.  Adaptive  Control  and  Guidance  for  Tactical  Missiles,  Volumes  I and  II. 

Reading,  Mass.,  30  June  1970.  (Report  no.  TR170-1.)  1970. 

®Naval  Weapons  Center.  Improved  AIM-7?  Guidance  law  GtiAy,  by  H.  D.  Nuffer.  China  Lake,  California,  NWC,  August 
1975.  (Working  papers.) 

^Air  Force  Systems  Command , Space  System  Division.  An  Optimum  Interception  T/xw  with  hounded  Control  in  Presence  of 
Noise,  bv  N.  E.  Nahl  and  D.  C.  Svorder.  AFSC,  Washington,  D.C.,  1967.  (Report  No.  SSD-TF-67-124.) 

10L.  A.  Stockum  and  F.  C.  Welmcr.  "Optimal  and  Suboptlmal  Guidance  for  a Short-Range  Homing  Hlssile,”  in  Inst. 
Elea.  Electron.  Eng.,  Trans.,  Aerospace  Electron.  Syst.,  Vol.  AES-12,  No.  3 (Hay  1976). 

11R.  G.  Cottrell.  "Optimal  Intercept  Guidance  for  Short-Range  Tactical  Missiles,"  AIAA  . Journal , Vol.  9,  No.  7 
(July  1971)  pp  1414-1415. 

*2R.  W.  Riche! . "Optimal  Torsi""'  Guidance  of  an  Alr-to-Surface  Missile,'  J.  Spacecraft,  Volume  5,  No.  6 
June  1968. 

12S.  M.  Bralnin  and  R.  B.  McGhee.  "Optimal  Biased  Proportional  Navigation,"  In  Inst.  Elec.  Electron.  Eng.,  Trans., 
Automat.  Cor.tr.,  Vol.  AC-13,  No.  4 (August  1968). 


- , j'i  <% 


tetaua-vm*?.  ■KaeLxcx-.R'  -t , y r^*r . (•.nr.AVE/.Sfe  t 


•'  -i«tf  ’«* 


KIKST?: 


NWC  TP  5913 


where 


A = commanded  missile  acceleration 
m 

c 

N = navigation  ratio 
a = missile/target  line  of  sight  rate 
aT  = target  acceleration 
• 

R = closing  velocity 


Such  a guidance  law  is  referred  to  in  Footnotes  1,  2,  7 and  8.  While 
the  derivation  of  Equation  1 can  be  mathematically  derived  in  a routine 
fashion  without  physical  insight,  it  is  also  intuitively  appealing. 

This  can  be  seen  by  studying  Figure  1 which  shows  a target  and  missile 
iri  the  horizontal  plane  in  two  different  positions  (Figure  la  and  b) . 

In  both  figures  the  target  begins  a hard  turn  at  time  t = tp  If  the 
missile  is  at  the  position  denoted  by  t = t^  in  Figure  la,  a change  of 
line  of  sight,  Ao^,  will  be  observed  in  a short  time  At.  On  the  other 
hand  if  the  target  is  already  at  the  position  denoted  by  t = tg  (Figure 
lb)  a much  larger  line-of-sight  angle,  A a2>  is  observed  in  a similar  time 
interval  At.  Noting  that  proportional  navigation  uses  a = as  its 

means  of  guidance,  it  is  seen  that  effectively  the  target  maneuver  is  not 
fully  detected  until  the  target  has  maneuvered  for  some  period  of  time. 

The  use  of  the  target  acceleration  term  a^  in  Equation  1 causes  the  missile 
to  make  a guidance  correction  as  soon  as  the  maneuver  starts.  As  the  simu- 
lation results  presented  later  in  this  report  indicate,  this  results  in  a 
more  direct  (shorter  length)  path  between  the  target  and  missile  and  a 
corresponding  shorter  time  of  missile  flight  than  that  obtained  using 
proportional  navigation. 

In  order  to  use  the  guidance  law  given  by  Equation  1,  it  is 
necessary  to  know  the  target  acceleration,  a,,,.  The  method  most  commonly 

X 

suggested  is  to  use  a Kalman  filter.  In  this  approach  the  target  is 
assumed  to  be  undergoing  continuous  random  maneuvers.  It  is  nec- 
essary to  have  an  estimate  of  the  target  maneuver  statistics  in  orde>. 
to  apply  this  technique.  In  addition,  use  of  the  basic  Kalman 
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filter  algorithm  assumes  a linear  model  of  the  missile/target  system 
which  implies  small  angle  approximations  in  the  kinematic  equations. 
Another  approach  to  estimating  the  target  acceleration  which  also  assumes 
small  angular  changes  in  the  missile/target  line  of  sight  is  to  relate 
the  target  acceleration  to  missile  acceleration,  line-of-sight  rate,  and 
acceleration  (see  footnote  7) 


a = Ra  - 2RCF  + a 
i m 


(2) 


where 


R = range 

a = missile  normal  acceleration 
m 

o = line-of-sight  acceleration 
• * 

While  R,  R,  a and  a^  can  be  readily  measured  by  existing  sensors  in 
missiles,  0 is  more  difficult  to  obtain.  Since  a can  be  expected  to  be 
noisy  in  the  actual  missile  hardware,  differentiation  of  a to  obtain  a 
is  undesirable.  A different  approach  to  estimating  target  acceleration 
is  presented  in  this  report.  It  is  not  necessary  to  make  small  angle 
approximations  or  to  have  knowledge  of  target  statistics  but  rather  that 
the  target  maneuver  can  be  approximated  by  an  arc  of  a circle  whose 
radius  is  changing  with  time. 

In  the  remaining  sections  of  this  report  the  target  velocity  and 
acceleration  estimator  is  derived  and  analyzed  after  implementing  it  in 
a missile  simulation.  A guidance  law  similar  to  that  given  by  Equation 
1 is  used  to  compare  trajectories,  time  histories,  and  miss  distances 
with  those  obtained  by  using  proportional  navigation. 
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DERIVATION  OF  TARGET  VELOCITY 
AND  ACCELERATION  ESTIMATOR 


In  this  section  of  the  report  the  basic  algorithm  for  estimating 
target  velocity  and  acceleration  is  presented.  Only  the  case  when  the 
target  and  missile  motion  is  constrained  to  a single  plane  is  considered 
here.  The  basic  assumptions  in  the  derivation  are  that  the  target  speed 
is  constant  and  ’ ■’t  the  target  is  either  not  maneuvering  or  is  turning 
in  a circular  arc.  Simulation  results  are  presented  later  in  this  report 
when  these  two  assumptions  are  not  met. 

Figure  2 shows  the  trajectories  of  the  missile  and  target  in  a 

common  plane  taken  to  be  the  horizontal  inertial  plane.  As  the  missile 

moves  from  the  position  denoted  by  (x  (t..),  y (t.))  to  (x  (t„),  y (t„)), 

ml  ml  m z m z 

the  target  moves  on  the  circle  from  the  position  denoted  by  0(t^)  to 
that  denoted  by  S^).  It  may  turn  out  that  the  target  is  not  maneuver- 
ing, i.e.,  the  trajectory  is  a straight  line  like  that  shown  between 
times  t = 0 and  t = t^  in  Figure  2.  However,  this  can  also  be  thought 
of  as  a circle;  the  radius  in  this  case  being  infinite. 


FIGURE  2.  Missile/Target  Kinematics  Used  in  Estimator  Derivation, 


* -j-y  ^ , 
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In  terms  of  the  missile  and  target  coordinates  (xm(t) , ym(t))  and 
(xT(t) , yT(t)),  respectively,  the  range  R(t)  between  target  and  missile 

is  given  by 


R(t)  = \ZVt}  + Ry(t) 


where 


\(t)  = xT(t)  - xm(t) 


Ry(t)  = yT(t)  - ym(t) 


Evaluating  Equation  3 at  times  t^  and  t2>  we  have  the  following 


expression 


R2(t2)  - R2(tL)  = R2(t2)  + R^(t2)  - RyCtp  - Rx(tx) 


Using  the  definition  of  a derivative,  we  also  have 


, , d R2 ( t)  _ lim  R2 ( t + h)  - R2(tj 

R(t)R(t)  jjt  2 h->o  2h 


In  Equation  6 let  t2  ® t + h and  t^  - t,  and  substitute  the  result  into 


Equation  7. 


li„  Ry<‘  + W - Ry<«  , lim  R*(t  + M ~ V£ 

R(t)R(t)  = 2h~  ' h+o  2h 


„ 4<‘> , d <(t) 

dt  2 dt  2 


Ry(t)Ry(t)  + Rx(t)Rx(t) 


i y»  *"**  **-  *\.’**HV  ' 


8 


or 


1 d_ 

1 + (R  /R  )2  dt 
y x7 


R2Cf  » R R - RR  (10) 

x y y x 

Equations  8 and  10  describe  the  .missile  and  target  regardless  of  the 
motion  of  either. 

At  this  point  the  assumption  of  circular  target  motion  and  constant 

target  speed  is  used.  As  indicated  in  Figure  3,  the  target  normal 

acceleration  a-r  is  directed  toward  the  center  of  the  circle  and  has 
n 

magnitude 


aT  - VT0  (11) 

n 

where 

7,j,  ■ magnitude  of  target  velocity 

0 = angular  rotation  rate  of  target 

a^  » target  acceleration  normal  to  its  velocity  vector 
n 

Integrating  Equation  11  between  t^  and  t2>  we  obtain 
aT  (t2  - tj)  - VT[0(t2)  - 9(t1)J 


or 


t1  + t2 


10 


(12) 
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where  t^  and  are  arbitrary  times  corresponding  to  different  positions 
of  the  target  along  the  circular  path. 


The  x and  y components  of  the  target  velocity  are  given  by 


xT(t)  = VT  sgn (g^)  sin  6(t) 


y (t)  = V sgn(0  ) cos  9(t) 


(13) 

(14) 


where 


xT(t) 


x component  of  target  velocity 
y component  of  target  velocity 
0(t)  = angular  position  of  target  (see  either  Figure  2 or  3) 


yT(t) 


(+i  if  gr  > o 

sgn (gT)  = 

(-1  if  9t  < o 

(15) 

it 

H 

aT  /9.807 

(16) 

n 


The  use  of  sgn(^)  is  to  keep  track  of  the  direction  of  rotation  about 


the  circle  shown  in  Figure  2.  Positive  acceleration  (g  > 0)  is  defined 


to  be  in  the  counterclockwise  direction  while  negative  acceleration 
(S’j  < 0)  is  defined  to  be  in  the  clockwise  direction.  From  Figure  3 we 


also  see  that 


R (t)  = R cos  a 
x 


(17) 


Ry(t)  = R sin  a 


(18) 
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Equations  8,  10,  13,  14,  17  and  18  provide  the  basis  for  the  target 
velocity  and  acceleration  estimator.  Equations  8,  10,  17,  and  18  can  be 
combined  to  give  the  following  equations. 


^-sin  a cos  a' 


cos  a sin  o 


1 

**•  *" 
R 

X 

*R a" 

m 

-1 

R 

L yJ 

-R  - 

(19) 


Inverting  and  solving  for  R and  R , we  have 

x y 


R = -Rcr  sin  a + R cos  a 
x 


R “ R0  cos  0 + R sin  0 

y 


(20) 

(21) 


By  definition 


• • • 


R = x_  - x 
x Jl  m 


(22) 


Ry  = yT  " ym 


(23) 


Rewriting  Equations  13  and  14  and  using  Equations  20  through  23,  we  have 
the  following: 


• • • ■ • 

XT  = s®n^T'  s*n  ® = xm  + R cos  0 ~ R0  sin  a 


(24) 


yT  " VT  sgn(^T)  cos  0 - ym  + R sin  0 + R0  cos  a 


(25) 


The  magnitude  of  the  target  velocity  is  then 
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vi  * \pr  * *i 


(26) 


v; 


. . • • .?  * 2 
= \/(xm  + R cos  a - Ra  sin  a)  . + (y  + R sin  o +Ra  cos  a) 


and  the  angular  position  (see  Figure  3)  of  the  target  is  given  by 


tan  6 = 


VT  sin  6 VT  sgr(0T)  sin  9 


VT  cos  0 VT  sgnCg^)  cos  6 


(27) 


Substituting  Equations  24  and  25  into  Equation  27  gives 


tan  9 


x + R cos  0 - Ra  sin  a 
m 


(28) 


y + r sin  a + Ra  cos  a 
ni 


or 


0(t)  = tan 


-1 


xm(t)  + R(t)  cos  a(t)  - R(t)a(t)  sin  a(t) 


ym(t)  + R(t)  sin  a(t)  + R(t)a(t)  cos  a(t) 


(29) 


Once  9(t)  is  determined  at  two  different  values,  t^  and  t^  of  time, 


they  can  be  substituted  along  with  VT  into  Equation  12  to  give  the 


target  acceleration  aj  normal  to  its  velocity  vector. 

n 


In  many  cases  it  is  not  necessary  to  compute  the  arc  tangent  as 
given  in  Equation  29;  an  approximation  works  equally  well.  If  we  let 


At  = - t^  * "sampling  interval" 


(30) 


we  can  use  Equa..ion  28  to  write  Equation  12  as 


aBKBfc  r isar.  r-m~ . & 
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aT 

n 


VT  0(t2)  - 9(tx) 
At 


_T  tan 
At 


-Yvv\  -i  /v^v 

Iwj  W 


aT  At 
n 


tan  tan 


-1  /VV^\  -1  /xT(tlA 

w/  why 


A trigonometric  identity  can  be  applied  to  the  right  side  of  Equation  32 
to  give  the  following 


aT  At? 
n 


tan  tan 


-1  /W' 


- tan  tan 


-1  /XT(V 


v*T<h> 


1 + tan 


-1  A 

l \W/J 


tan  tan 


XT^2^ 


yT(t2)  yT(tj) 


, ww 

yT(t2)  y^'tj) 


If  we  assume  that  we  can  replace  tan  l-y 


following 


^ At 

T 

n 


aT  At 
n 


, we  have  the 


f b(t2)  - b^)  J 
(1  + b(t2)  b(t1)j 


M.  „ _ 


14 


}< 


NWC  TP  5913 


where  . 


xT(t) 

b(t)  = — 

yT(t) 


The  approximation  that 


(35) 


(3o) 


though  not  necessary,  would  seem  to  be  realistic  for  most  applications. 

For  example,  if  the  target  is  pulling  10  ^'s,  (a^  a 98.07  meters  per 

n 

2 

second  squared  (m/s  )),  the  target  velocity  is  200  m/s  (656  fps)  and 
the  acceleration  is  to  be  estimated  every  tenth  of  a second  (At  = 

0.1  sec),  we  have 


aT  At 

— = 0.049035  radians 

T 


while 


tan 


1=  0.049074  radians 


As  Equation  36  indicates,  the  smaller  the  target  acceleration,  the 
shorter  the  time  interval  between  estimation  samples  and  the  higher  the 
target  velocity,  the  better  the  approximation  in  Equation  34  will  be. 
All  of  the  simulation  results  presented  later  in  this  report  estimate 
the  target  acceleration  by  using  Equation  34. 
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While  Equations  31  or  34  can  be  used  to  estimate  the  target 
acceleration  normal  to  the  target  velocity  vector,  it  is  not  this 
acceleration  that  is  used  in  the  guidance  law  given  by  Equation  1 but 
rather  the  acceleration  perpendicular  to  the  lint  of  sight  between 
target  and  missile.  This  component  of  the  normal  acceleration  can  be 
determined  from  the  geometry  shown  in  Figure  4.  The  component  of  target 
acceleration  perpendicular  to  the  line  of  sight  is 


sin  (o  + 0) 


(37) 


a_  can  thus  be  determined  by  using  Equation  29  to  determine  9 and 
TL 

either  Equation  3x  or  34  to  determine  aT  while  the  line-of-sight  angle 

n 

0 is  assumed  to  be  measurable.  The  purpose  of  the  absolute  value  of  aT 

n 

instead  of  just  a^  is  to  account  for  the  direction  of  rotation  of  the 
n 

target  and  to  assign  a direction  to  a . 

U 


FIGURE  4.  Geometry  for  Determining  Target  Acceleration  a 
Perpendicular  to  Line  of  Sight. 


16 


NWC  TF  5913 

An  alternate  expression  for  aT  can  be  obtained  by  expanding 

1 

Equation  37: 

a = |a  | (sin  0 cos  9 + sin  0 cos  0)  (38) 

J.  < X 

X n 

Substituting  Equations  13,  14,  17,  and  18  into  Equation  38,  gives 


ri. ' |aiJ  ( 


R y_  x,,,  R 

y JT  , T x 

R VT  sgn(0T)  VT  sgn(g'T)  R 


|aT  | sgn(aT  )(R  yT  + x^) 
1 n1  n y 


V 


aTn(RyyT  + RxXT) 

V 


(39) 


In  this  form  a-p^  can  be  obtained  without  computing  the  arc  tangent 
shown  Jn  Equation  29  as  long  as  the  approximation  in  Equation  36  is 
accepted.  The  components  of  target  velocity,  xT»  y^,  shown  in 

Equation  39  are  given  by  Equations  24  and  25. 


SUMMARY  OF  ESTIMATOR  ALGORITHM 

The  method  of  determining  target  velocity  and  acceleration  for  a 
target  following  a circular  path,  as  derived  ir  the  previous  section, 
provides  the  basis  of  the  general  target  velocity  and  acceleration 
estimator.  Using  those  results,  the  algorithm  shown  in  Figure  5 is 
obtained.  In  this  figure  the  symbol  ~ is  placed  over  the  estimated 
variable,?  to  distinguish  them  from  the  true  variables. 
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FIGURE  5.  Target  Velocity  and  Acceleration  Estimator  Algorithm. 
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One  of  che  inputs  to  the  estimator  algorithm  is  the  time  interval 
At  between  estimates  of  the  target  acceleration.  From  the  point  of 
accuracy  it  would  be  desirable  for  this  interval  to  be  small  while  from 
a computational  burden  point  it  would  be  desirable  for  At  to  be  large. 
Simulation  results,  some  of  which  are  presented  later  in  this  report, 
indicate  that  a value  of  At  =■  0.1  sec  provides  reasonable  accuracy  in 
the  estimate.  Additional  inputs  to  the  estimator  include  range,  range 
rate  (closing  velocity) , line-of-sight  angle  and  rate  as  well  as  components 
of  missile  velocity.  Note  that  the  line-of-sight  angle  and  components 
of  missile  velocity  are  referenced  to  an  inertial  system  and  not  to  the 
body  of  the  missile.  The  components  of  target  velocity,  x_,  y_,  are 
then  estimated  using  these  inputs  and  Equations  24  and  25.  The  magnitude 
VT  of  the  target  velocity  is  then  obtained  from  Equation  26. 

In  order  to  estimate  the  target  acceleration  it  is  necessary  to 
estimate  and  store  the  components  of  estimated  target  velocity  at  two 
times  separated  by  the  sampling  interval  At.  Equations  34  and  35  are 
then  used  with  t^  * t - At  and  t2  = t to  estimate  the  target  accelera- 

ation  normal  to  its  velocity  vector.  Finally,  Equation  39  is  used  to 
estimate  the  target  acceleration  normal  to  the  line  of  sight  between 
target  and  missile. 


EVALUATION  OF  ESTIMATOR 

In  order  to  gain  some  insight  into  how  well  the  estimator  works  and 
what  can  be  gained  in  the  way  of  improved  missile  performance  by  using 
such  an  estimator,  a simple  model  of  a missile  and  target  was  constructed 

and  simulated  on  a digital  computer.  The  accuracy  of  the  estimator  was 
observed  by  comparing  estimated  values  of  target  velocity  and  accelera- 
tion with  exact  values.  This  was  done  under  conditions  not  satisfying 
the  assumptions  under  which  the  estimator  was  derived  as  well  as  under 
those  conditions  for  which  it  was  derived.  Time  histories  of  missile 
acceleration  and  missile/target  trajectories  were  compared.  The  variable 
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in  the  comparison  was  the  guidance  law,  i.e. , missile  performance  was  com- 
pared with  and  vithout  the  target  acceleration  estimator  in  the  guidance 
law.  The  two  guidance  laws  compared  were 


m 


NRcr 


(proportional  navigation) 


(40) 


and 


A «'N  Ro  - ,S(t  )a_ 
m_  2 v go  T 


(41) 


where 


A “ commanded  missile  acceleration  measured  in  the  wind 
c axis  system 

N * navigation  ratio  * 3 

6 ■ missile/target  line-of-sight  rate 

a = estimated  target  acceleration  perpendicular  to  the 
missile/target  line  of  sight 

R * closing  velocity 


S<V 


1 if  t > t sec 

90  - fl^min 
0 if  t < t sec 
90  go  mln 


t * R/R  (approximate  "time  to  go"  until  intercept) 

g° 


(42) 

(43) 


This  is  essentially  a comparison  of  proportional  navigation  with  pro- 
portional navigation  biased  by  the  estimated  target  acceleration.  A 
comparison  of  Equations  1 and  41  reveals  the  addition  of  a switch  S(t  ) 
in  Equation  41  which  switches  the  estimator  out  of  the  guidance  law  when 
the  approximate  "time-to-go"  until  intercept,  t , is  less  than  some 

small  value  t_  . The  use  of  this  switch  is  included  because  in 
Strain 

some  cases  the  estimate  of  target  acceleration  became  poor  when  the  time 
to  go  until  intercept  became  less  than  0.5  sec  (see  Figure  14  for 
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example) . This  decrease  in  accuracy  of  the  estimator  is  thought  to  be 
due  to  the  inherent  instability  of  the  kinematic  equations  as  the  range 
R goes  to  zero.  In  the  simulation  examples  presented  later,  values  of 
t^o  of  1 and  2 sec  were  used.  In  all  example  trajectories  presented 

below,  with  the  exception  of  Example  5,  the  target  acceleration  estimate 

/v 

a used  in  Equation  41  was  updated  every  0.1  sec.  Between  updates 
rX 

the  target  acceleration  estimate  is  held  constant  at  the  value  of  the 
most  recent  estimate.  The  constant  1/2  appearing  in  Equation  41  results 
from  theoretical  guidance  law  studies  (see  Footnotes  1,  7,  and  8)  based 
on  simplifying  assumptions.  No  attempt  was  made  to  optimize  this 
parameter.  Its  value  could  be  determined  by  a Monte-Carlo  study. 

SIMULATION  MODEL 

The  missile  and  target  were  modeled  as  point  masses.  Only  motion 
in  a horizontal  plane  was  considered.  A constant  missile  velocity  of  1 
kilometer  per  second  (.54  nautical  mile  per  second)  was  assumed. 

Target  motion  varied  from  one  missile/target  encounter  to  the  next  as 
described  in  the  example  computer  simulations  described  below.  The 

commanded  missile  accelerations  as  given  by  Equations  40  or  41  were  i 

} 

applied  to  an  acceleration  command  limiter  and  second  order  airframe 
model  as  shown  in  Figure  6.  Numerical  values  of  the  parameters  appear- 
ing in  Figure  6 correspond  to  a missile  currently  under  investigation 
and  are  discussed  further  in  the  examples  shown  below.  Missile  accel- 
eration A is  converted  from  wind  axes14  to  inertial  coordinates  and 
m 

integrated  to  give  velocity.  Since  the  magnitude  of  the  missile  velocity 
will  tend  to  increase  slightly  due  to  numerical  integration  errors,  the 
magnitude  of  the  velocity  is  normalized  to  the  assumed  constant  value  of 
1 kilometer  per  second  (.54  nautical  mile  per  second).  The  resultant  com- 
ponents of  missile  velocity  are  then  integrated  to  provide  missile  position 
so  that  the  relative  position  of  the  missile  and  target  can  be  computed. 

14  1.  Kevin.  Dyrmiet  ef  Flight.  John  Wiley  end  done,  We*  York,  1959. 
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ACCELERATION  LIMITER 


AUTOPILOT-SERVO-AIRFRAME 


U>n  - AUTOPILOT-SERVO-AIRFRAME  NATURAL  FREQ  (RAD/SEC) 
f « AUTOPILOT-SERVO-AIRFRAME  DAMPING  RATIO  » 0.707 
S - LaPLACE  TRANSFORM  VARIABLE 

FIGURE  6.  Missile  Pitch  Plane  Model. 


It  is  noted  that  the  use  of  the  target  acceleration  estimator  given 
in  this  report  is  not  restricted  to  or  derived  on  the  basis  of  this 
simplified  model,  but  rather  it  is  being  used  as  a simple  tcol  to  study 
its  performance.  On  the  following  pages  a number  of  example  missile/ 
target  encounters  which  compare  missile  performance  with  and  without  the 
target  acceleration  estimator  are  simulated  and  discussed. 

Example  1.  In  this  example  the  missile  and  target  are  initially 
separated  by  22  kilometers  (11.9  nautical  miles)  in  a head-on  situation. 
Time  t * 0 sec  is  taken  to  be  the  time  of  the  beginning  of  the  endgame 
encounter.  Numerical  values  of  the  missile  parameters  shown  in 
Figure  6 are  given  in  Table  1.  As  shown  in  Figure  7,  for  6 sec  after 
the  beginning  of  the  terminal  encounter  the  target  follows  a straight 
line  trajectory  and  then  begins  a 5-g  turn  in  the  counterclockwise 
direction.  Figure  7 compares  the  trajectories  of  the  missile  with  the 
use  of  the  target  estimator  and  without  it  (proportional  navigation) . 

As  can  be  seen  from  this  figure,  the  use  of  the  target  acceleration 
estimator  provides  a more  direct  or  shorter  trajectory  to  the  target  due 
to  its  ability  to  compensate  for  the  maneuver.  Both  missiles  guided 
successfully  to  the  target,  however.  Figure  8 compares  exact  missile 
velocity  and  acceleration  with  those  computed  using  the  estimation 
technique  given  in  this  report.  The  quality  of  the  estimation  is  seen 
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TABLE  1.  Missile  and  Target  Simulation  Parameters. 

Example 


\ Parameters 

1 

. 2 

3 

4 



5 

| 

Missile 

t 

l Airframe  natural  frequency 

’ to  . seo“> 

t n 

6.4 

1.28 

6.4 

6.4 

6.4 

Airframe  damping  ratio,  5 

0.707 

0.707 

0.707 

0.707 

0.707 

Velocity,  meter /second 

1000.0 

1000.0 

1000.0 

1000.0 

1000.0 

Altitude,  kilometer 

6.0 

23.0 

6.0 

6.0 

6.0 

t Maximum  acceleration 

15.0 

5.3 

15.0 

15.0 

15.0 

{ (Amax  in  Pi£ure  6) 

, Estimator  update 

0.1 

0.1 

0.1 

0.1 

0.01 

; interval  At,  second 

f Target 

Velocity,  meter /second 

379.0 

i036.0 

379.0 

Sinusoidal 

Sinusoidal 

} 

(Mach  1.2) 

(Mach  3.5) 

(Mach  1.2) 

Altitude,  kilometer 

6.0 

23.0 

6.0 

6.0 

6.0 

; Maneuver  magnitude,  g1 a 

5.0 

2.0 

5.0 

Variable 

Variable 

i Maneuver  initiation, 

6.0 

Variable 

6.0 

0 

0 

second 

| Maneuver  time  constant, 

0 

0 

1.0 

• 4 • 

! • second 

; Range  from  missile  at 

22.0 

28.0 

22.0 

22.0 

22.0 

f launch,  kilometer 



10.0  r 


7.5  U 


a.  5.0 

ui 

I- 

UJ 

2 

o 

d 2.5 
X 

w* 

X 

5 0.0 


-2.5 


-5  0 


ALTITUDE  - 6 KILOMETERS 
MISSILE  VELOCITY  « 1000.0  m/s 
TARGET  VELOCITY  - 379.0  m/s  (MACH  1.2) 

TARGET  MANEUVER  « 5-g  TURN  INTO  MtSSILE  STARTING  AT  t»  6.0  SEC 
INITIAL  RANGE  - 22  KILOMETERS  (11.9  NAUTICAL  MILES) 


WITHOUT  TARGET  ACCELERATION 


t-6  SEC  t *=  0 


MISSILE 

MOTION 


t = 6 SEC  

*'  ^r.' target 

WITH  TARGET  ACCELERATION 


MOTION 


0.0  2.5  5.0  7.5  10.0  12.5 

X AXIS,  KILOMETERS 


15.0 


17.5 


20.0 


22.5 


FIGURE  7.  Missile  and  Target  Trajectory  in  the  Horizontal  Plane 
for  Example  1, 


23 


NWC  TP  5913 


ALTITUDE  = 6 KILOMETERS 
TARGET  MANEUVER  - -5b 


-ACTUAL  ACCELERATION 


•COMMANDED  ACCELERATION 


(/> 

o>  10.0 


, TARGET  ACCELERATION  ESTIMATOR 
SWITCHED  OUT  OF  GUIDANCE  LAW 


ui  0.0 

O 


< 

S 

cc 

§1  10.0 


TARGET  MANEUVER 
INITIATION 


TIME  UNTIL  INTERCEPT,  SEC 
(a)  WITH  TARGET  ACCELERATION 


ALTITUDE  = 6 KILOMETERS 
TARGET  MANEUVER  - -50 


— ACTUAL  ACCELERATION  *=  A. 


• COMMANDED  ACCELERATION  » A_ 


TARGET  MANEUVER 
INITIATION 


10.0  12.6 


15.0  17.5 


TIME  UNTIL  INTERCEPT,  SEC 

(b)  WITHOUT  TARGET  ACCELERATION  (PROPORTIONAL  NAVIGATION) 


FIGURE  9.  Missile  Acceleration  Profiles  for  Example  1. 
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to  be  very  good.  Figure  9 compares  missile  acceleration  with  and  without 
the  target  acceleration  estimator.  The  acceleration  A and  Am  shown 
in  this  figure  refer  to  those  shown  in  Figure  6.  As  shown  in  Figure  9, 
the  target  begins  to  maneuver  when  the  time  until  intercept  is  approxi- 
mately 11  sec. 

Comparing  missile  acceleration  profiles  with  and  without  the  target 
acceleration  estimator  show  that  the  missile  with  the  target  acceleration 
estimator  effectively  reverses  its  direction  approximately  2 sec  sooner 
than  the  one  without  the  estimator.  This  ability  to  detect  the  maneuver 
earlier  results  in  its  more  direct  trajectory  for  the  missile  with  target 
acceleration  estimator  as  shown  in  Figure  7,  Again  referring  to  Figure  9, 
an  abrupt  change  in  the  commanded  acceleration  occurs  in  the  missile  em- 
ploying the  target  estimator  when  the  time  until  intercept  is  approxima- 
tely 2 sec.  This  is  due  to  the  estimator  being  switched  out  of  the  guid- 
ance law  given  by  Equation  41  when  T = 2 sec  in  Equation  42. 

* min 

This  example  was  repeated  using  a maximum  missile  lateral  accelera- 
tion of  5 g instead  of  15  g.  The  missile  using  proportional  navigation 
missed  the  target  by  440  meters  while  the  missile  employing  the 
estimated  target  acceleration  missed  by  less  than  4 meters. 

Example  2.  The  ability  of  the  missile  with  the  target  acceleration 
estimator  to  more  quickly  change  its  course  when  a maneuver  starts  and 
to  follow  a more  direct  trajectory  as  indicated  in  Example  1 was  thought 
to  be  of  particular  advantage  at  high  altitude  where  missile  maneuverability 

is  low  and  airframe  response  is  slow.  In  this  example,  the  altitude  of 
missile  and  target  is  changed  to  23  kilometers.  Maximum  missile 
acceleration  is  limited  to  5.3  g' s while  the  target  maneuver  is  taken  to 
be  2 g'  & as  indicated  in  Table  1.  Figure  10  shows  the  trajectories  of 
target  and  missiles  with  and  without  the  target  acceleration  estimator 
for  t _>  6 sec  where  t = 6 sec  is  the  time  at  which  the  maneuver  starts. 

As  in  Example  1 the  use  of  target  acceleration  estimation  provides  a 
more  direct  path  to  the  target.  In  this  case  the  missile  without  the 
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target  acceleration  estimator  missed  the  target  by  60  meters  while  the 
missile  with  the  estimator  came  within  5 meters  of  the  target.  (In  all 
examples  considered  in  this  report  the  simulation  was  stopped  whenever 
the  missile  came  within  5 meters  of  the  target.)  The  corresponding 
missile  acceleration  time  histories  are  shown  in  Figure  11.  From  this 
figure  it  can  be  seen  that  because  the  missile  without  the  acceleration 
estimator  is  acceleration  limited  the  last  4 sec  of  flight,  it  is  unable 
to  reach  the  target.  However,  the  missile  with  the  target  acceleration 


estimator  reverses  its  acceleration  sooner  after  the  maneuver  starts  and 


is  able  to  guide  to  the  target.  Figure  12  shows  the  comparison  between 
true  and  estimated  target  acceleration  and  velocity;  again  the  agreement 
is  very  good. 


ALTITUDE  - 23  KILOMETERS 

MISSILE  VELOCITY  - 1000.0  m/s 

TARGET  VELOCITY  = 1036.0  m/s  (MACH  3.5) 

TARGET  MANEUVER  - 2-g  TURN  INTO  MISSILE  STARTING 
AT  t - 6 SEC 

INITIAL  RANGE  - 28  KILOMETERS  (15.1  NAUTICAL  MILES) 


| 2.0 
O 


TARGET 

MOTION 


52  0.0 

X 

< 


MISSILE 

MOTION 
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FIGURE  10.  Missile  and  Target  Trajectories  in  the  Horizontal  Plane 
for  Example  2. 
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The  previous  paragraph  was  concerned  with  a target  maneuver  start- 
ing 6 sec  from  the  beginning  of  the  endgame.  In  Figure  13  are  shown 
miss  distances  for  target  maneuvers  of  the  same  acceleration  magnitude 
but  starting  at  different  times.  As  this  figure  indicates,  significant 
miss  distance  reduction  is  possible.  In  the  case  when  the  maneuver 
starts  12  sec  after  the  terminal  encounter  begins,  the  miss  distances 
with  and  without  the  target  acceleration  estimator  are  the  same.  This 
is  due  to  the  maneuver  being  started  when  the  time  until  intercept  is 
less  than  2 sec.  In  this  example,  2 sec  until  intercept  was  also  chosen 
to  be  the  time  at  which  the  estimator  was  switched  out  of  the  guidance 
law  given  by  Equation  41.  Thus,  for  the  case  when  the  maneuver  starts 

at  12  sec  into  the  endgame,  the  target  acceleration  estimator  was 
actually  not  used  and  the  missile  response  is  the  same  as  if  propor- 
tional navigation  had  been  used,  leading  to  the  same  miss  distance. 

Example  3.  In  the  previous  examples  it  was  assumed  that  the  target 
was  moving  with  constant  velocity  and  instantaneously  changed  it  re- 
acceleration from  one  constant  level  to  another  (see  Figure  12  for 
example).  These  two  assumptions  ar*  the  basic  assumptions  under  which 
the  estimator  was  derived,  and  we  would  expect  good  results  under  these 
conditions.  In  this  example  the  second  assumption  is  removed,  but  the 
first  assumption  (constant  target  velocity)  is  retained.  More  speci- 
fically, the  target  is  switched  from  0 to  -5  ^'s  starting  6 sec  after 
the  endgame  begins  but  changes  level  exponentially  (see  Figure  14). 
target  acceleration  is  modeled  by 

/ -(t-6)/T  \ 

gT  = gT  \1  - e ) for  t > 6 sec  (44) 

c 

= target  acceleration  (gravity  units) 

= commanded  target  acceleration  = -5  g's 

= target  time  constant  * 1 sec 


Thus  the 


where 
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FIGURE  14.  Comparison  Between  Estimated  and  Actual  Target 
Acceleration  und  Velocity  for  Example  3. 
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Simulation  parameters  are  given  in  Table  1.  This  example  differs 
from  Example  1 only  in  that  the  target  time  constant  tt  is  1 sec  in  this 
example  and  0 sec  in  Example  1.  Trajectories  of  the  target  and  the 
missile  with  and  without  the  target  acceleration  estimator  are  shown  in 
Figure  15.  Again,  the  missile  that  employs  the  acceleration  estimator 
follows  a more  direct  (shorter)  trajectory  to  the  target  than  the  same 
missile  without  the  estimator. 


FIGURE  15.  Missile  and  Target  Trajectories  in  the 
Horizontal  Plane  for  Example  3. 


The  accuracy  of  the  estimator  for  this  example  is  indicated  in 
Figure  14.  The  accuracy  of  the  target  velocity  is  nearly  perfect.  This 
follows  since  the  target  has  constant  velocity  and  this  is  one  of  the 
two  assumptions  in  the  derivation  of  the  estimator.  Also,  the  estimator 
determines  the  target  velocity  first  and  then  the  target  acceleration 
(not  vice  versa) , so  even  though  the  target  acceleration  is  not  of  the 
form  that  was  used  to  derive  the  estimator  the  target  velocity  estima- 
tion is  not  affected.  The  target  acceleration  is  still  good  although  it 
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tends  to  lag  in  time  the  true  value  at  the  time  the  maneuver  begins  (see 
Figure  14  with  "time  until  intercept"  ■ 11  sec).  Also,  the  estimate 
deteriorates  the  last  fraction  of  a second.  Since  the  estimator  was 
switched  out  of  the  guidance  law  (Equation  41)  when  the  time  until 
intercept  was  less  than  1 sec,  this  had  no  effect  on  the  missile  guidance. 

Example  4.  In  this  example  the  target  is  modeled  in  such  a way 
that  neither  of  the  assumptions  used  in  the  derivation  of  the  estimator 
is  satisfied.  The  coordinates  xT(t),  yT(t)  of  the  target  are  taken  to 
be 


*T(t)  • i:T(0)  - V t 

(45) 

yT(t>  ’ ymax  81,1  V 

(46) 

where 

V * magnitude  of  x component  of  target  velocity  * 379  meter/sec 
XT 

y = maximum  y component  of  target  position  =*  2000.0  meters 

tn&x 

» frequency  of  oscillation  * 0.24257  rad/sec 

Substituting  Equation  45  into  46  gives 

yT(t)  - ymax  8i"j“(*T<0>  " *T(t>)] 

(47) 

where 


“ ’ VV*T 


(48) 


0.000640  rad/meter 


(49) 
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Thus,  by  Equation  47  the  target  follows  a sinusoidal  trajectory  in  the 
horizontal  x,y  plane. 

Differentiation  of  Equations  45  and  46  gives 

xT(t)  - -Vx  (5C 


Ve)  ' W%  <=°»[“(*t<0)  " xi<c>)] 


V * Fl *2 
T WxT+yT 


\ \A  + coa2[<o(3cT(0)  - xT(t))J 


where 


H 


m 


VT  ■ target  velocity 


Thus  the  target  velocity  varies  between 


V * 379  meter/ sec 
XT 


Vx  \ fl  + y 2 m2  “ 615*6  meters/sec 
T V ymax 

Using  a result  from  Footnote  15,  the  target  acceleration  normal  to  its 
path  is  given  by 

yTxT  - xTyT 


P.  I.  Cook  and  D.  N,  Pierca.  Itynatioe,  Intarnational  Tax t book  Co«pai»y,  1960,  pagt  13. 
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Differentiation  of  Equations  50  and  51  and  substitution  into  Equation  53 
gives 


ymaxw2vxT  sln[W(x(0)  - xT(t))] 
Fn  JT+  yJxU)2  cos2[m(x(0)  - xT(f| 


(54) 


Thus  the  target  normal  acceleration  varies  between 


2 2 2 
±y  to  V = ±117.7  meters/sec  = ±12  <7*s 
'max  x^  a 


From  Equations  52  and  54  it  can  be  seen  that  the  assumptions  used  to 
derive  the  target  velocity  estimator,  i.e. , constant  velocity  and 
piecewise  constant  normal  acceleration,  are  not  satisfied. 

Figure  16  compares  the  trajectories  of  target  and  missile;  the 
latter  with  and  without  the  target  acceleration  estimator.  As  in  the 
other  examples,  use  of  the  target  acceleration  estimator  provides  a 
more  direct  trajectory  to  the  target.  Figure  17  compares  estimated  and 
actual  target  acceleration  and  velocity.  While  the  target  acceleration 
estimate  is  certainly  not  perfect,  the  estimated  acceleration  provides 
enough  bias  in  the  guidance  given  by  Equation  41  to  reduce  the  miss 
distance  from  30  meters  using  proportional  navigation  (i.e.,  no  target 
acceleration  estimate)  t;o  8 meters  with  the  target  acceleration  estimate. 
Figure  18  compares  missile  acceleration  with  and  without  the  target 
acceleration  estimator.  As  has  been  observed  in  many  simulation  runs, 
there  is  a tendency  for  the  missile  using  the  target  acceleration  esti- 
mator to  be  less  g- limited  (i.e.,  require  less  available  acceleration 
normal  to  its  flight  path)  than  the  missile  using  proportional  naviga- 
tion. A discontinuity  in  the  commanded  missile  acceleration  occurs  at 
approximately  1.25  sec  until  intercept  (see  Figure  18).  This  is  due 
to  the  target  acceleration  estimator  being  switched  out  of  the  guidance 
law  given  by  Equation  41  when  the  estimated  time  until  intercept  (as 
given  by  Equation  43)  is  less  than  1 sec. 
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FIGURE  16.  Missile  and  Target  Trajectories  in  the  Horizontal  Plane 
for  Example  4. 


Example  5.  This  example  differs  in  only  one  way  from 
Example  4.  In  Example  4 (as  well  as  all  the  other  previous  examples) 
the  time  interval  between  which  samples  of  range,  closing  velocity, 
line-of-sight  angle,  line-of-sight  rate,  and  missile  velocity  which  are 
used  to  estimate  the  target  acceleration  was  taken  to  be  0.1  sec.  In 
this  example  the  time  interval  is  decreased  by  a factor  of  10  to  0.01  sec 
to  see  if  the  estimated  acceleration  as  given  by  Figure  17  would  be 
improved  by  increasing  the  sampling  rate.  Figure  19  shows  the  comparison 
between  estimated  and  true  target  acceleration  and  velicicy  for  the 
smaller  sampling  interval  (At  = 0.01  sec).  Comparison  of  these  rf suits 
with  those  given  in  Figure  17  which  uses  the  longer  sample  interval 
(At  =0.1  sec)  reveals  little  improvement.  Consequently,  the  missile 
trajectory  and  acceleration  time  history  for  this  example  are  essentially 
the  same  as  those  shown  in  Example  4. 
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FIGURE  17.  Comparison  Between  Estimated  and  Actual  Target 
Acceleration  and  Velocity  for  Example  4 (Sampling  interval 
At  * 0.10  sec) . 
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O.L  2.5  5.0  7.5  10.0  12.5  15.0  17.5 
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(a)  WITH  TARGET  ACCELERATION. 


(b)  WITHOUT  TARGET  ACCELERATION  (PROPORTIONAL  NAVIGATION). 
FIGURE  18.  Missile  Acceleration  Profiles  for  Example  4. 


| ] FIGURE  19.  Comparison  Between  Estimated  and  Actual  Target 
B * Acceleration  and  Velocity  for  Example  5 (Sampling  interval 
Kt  | At  * 0.01  sec). 
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SUMMARY  AND  CONCLUSIONS 


In  this  report  a simple  algorithm  is  presented  for  estimating  target 
velocity  and  acceleration.  The  algorithm  uses  as  inputs  the  measure- 
ment of  range,  range  rate,  line-of-sight  angle,  line-of-sight  rate,  and 
missile  velocity.  The  estimated  target  acceleration  is  then  included  in 
a missile  guidance  law  and  the  effectiveness  of  such  a missile  compared 
to  a similar  one  using  proportional  navigation. 

The  primary  advantage  of  a missile  using  a target  estimator  is  that 
it  provides  a more  direct  (shorter)  trajectory  between  missile  and 
target  due  to  its  ability  to  correct  for  target  maneuvers  more  quickly. 
When  proportional  navigation  is  employed,  there  is  a tendency  for  the 
missile  to  be  acceleration  limited  (i.e. , pulling  maximum  g' s)  near  the 
end  of  its  trajectory  as  the  missile  tries  to  compensate  for  the 
maneuver.  The  use  of  the  target  acceleration  estimator  in  the  guidance 
law  provides  a means  of  correcting  for  the  target  maneuver  more  quickly 
and  thus  fewer  target  misses  due  to  missile  gr-limiting.  This  would  seem 
to  be  of  particular  advantage  at  high  altitudes  where  missile  maneuver- 
ability is  limited  and  airframe' response  is  slow  (see  Figure  13).  An 
additional  advantage  of  the  use  of  target  acceleration  estimation  in  the 
guidance  law  is  a slightly  shorter  time  of  flight  (0.2  to  0.5  sec  for  the 
examples  considered  in  this  report)  until  intercept  with  the  target. 

The  price  to  be  paid  for  target  acceleration  estimation  is  in  terms  of 
additional  sensors  and  system  complexity.  A comparison  of  Equations 
40,  41,  and  the  target  acceleration  estimator  as  given  by  Figure  5 
reveals  that  in  order  to  estimate  target  acceleration  the  additional 
variables  (i.e.,  those  variables  in  addition  to  closing  velocity  and 
line-of-sight  rate  used  in  proportional  navigation)  that  need  to  be 
measured  are  range,  line-of-sight  angle,  and  missile  velocity.  In 
addition,  a small  computer  is  necessary  in  order  to  implement  the 
acceleration  estimator  as  diagrammed  in  Figure  5. 
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The  estimation  algorithm  presented  here  appears  to  provide  a fairly 
good  estimate  of  target  acceleration  for  the  noise-free  cases  considered. 
In  the  actual  missile  one  can  expect  incorrect  measurements  of  range, 
range  rate,  level-of-sight  rate,  etc.,  needed  to  implement  the  algorithm 
presented  in  this  report.  These  measurement  errors  can  result  from 
glint,  clutter,  receiver  noise,  eclipsing,  countermeasures,  etc.  The 
sensitivity  of  the  estimator  to  these  noise  sources  was  beyond  the  scope 
of  this  study. 

Only  one  guidance  law  (i.e.,  the  one  given  by  Equation  41)  was 
investigated  in  this  study.  Also,  the  estimator  given  here  appears  to 
estimate  target  velocity  more  accurately  than  target  acceleration  (see 
Figure  17  for  example) . Unfortunately,  there  seems  to  be  little  in  the 
way  of  guidance  laws  based  on  target  velocity  rather  than  on  target 
acceleration.  As  mentioned  previously,  che  weighting  factor  of  1/2 
appearing  in  Equation  41  results  from  optimal  control  theory  based  on 
simplifying  assumptions.  It  is  not  known  how  sensitive  terminal  dis- 
tance is  to  this  parameter.  The  method  of  switching  the  target  accelera- 
tion estimator  out  of  the  guidance  law  (Equations  41  and  42)  produces  an 
abrupt  change  in  the  commanded  acceleration  (see  Figure  18  for  example). 
It  may  be  more  desirable  to  gradually  switch  the  acceleration  out  of  the 
guidance  law  than  to  use  Equation  42.  In  an  analog  circuit  a "soft" 
switch  can  be  implemented  using  a low-pass  filter  in  conjunction  with 
Equation  42.  Digitally,  it  can  be  done  by  replacing  Equation  42  with 


) = 


1 


isL 


+ rt 


n 

9° 


where  n and  r are  positive  numbers  which  can  be  chosen  to  provide  a less 
abrupt  switch  than  in  Equation  42. 

Lastly,  the  estimator  given  in  this  report  was  derived  under  the 
assumption  of  target  and  missile  coplanar  throughout  the  endgame.  While 
it  is  felt  that  the  generalization  from  two  to  three  dimensions  would  be 
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straightforward,  this  has  not  been  done  at  this  time  but  would  be 
necessary  if  such  an  estimator  were  to  be  actually  implemented  in  a 
missile. 

While  this  report  is  mainly  concerned  with  using  a target  velocity 
and  acceleration  estimator  in  a missile,  another  possible  application 
would  be  in  the  launch  aircraft  itself  as  a missile  fire  control  indi- 
cator. A proposed  method  for  indicating  when  the  missile  should  be 
launched  uses  precomputed  launch  regions  (launch  envelopes)  that  indi- 
cate to  the  pilot  that  the  geometry  relating  target  and  missile  is  such 
that  the  missile  has  a good  chance  of  intercepting  the  target.  Whether 
or  not  a missile  can  successfully  interrupt  a target  depends  to  a large 
degree  on  target  velocity  and  acceleration.  The  use  of  a target  acceler- 
ation and  velocity  estimator  such  as  the  one  proposed  in  this  report 
could  be  used  to  generalize  the  launch  aircraft  fire  controller  to  in- 
clude the  case  of  maneuvering  targets. 
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